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ABSTRACT 


A  paLtern  ayntnesis  method  is  presented  that  has 
approximative  characteristics  similar  to  that  of  a  Fourier 
aeries  but  is  computationally  much  simpler,  and  conse¬ 
quently  results  in  a  higher  over-all  accuracy.  Error 
estimates  and  convergence  measures  are  presented,  which 
compare  the  above  synthesis  to  the  Fourier  approximation 
and  the  Woodward  method. 


iii 


I.  INTRODUCTION 

The  two  moat  important  requirements  in  pattern  synthesis  are 
minimization  of  error  and  the  ease  of  computation  In  obtaining  the  ex¬ 
citation  coefficients.  Sometimes  one  of  the  two  has  to  be  sacrificed 
for  the  other.  Thus,  of  the  two  most  prevalent  methods,  the  Fourier 
expansion  is  the  synthesis  procedure  that  results  in  a  known  error  which 
is  also  minimized  in  the  mean  square  sense.  The  difficulty  in  obtaining 
the  excitation  coefficients  is  often  quite  serious,  particularly  if  one  has 
to  include  the  element  factor  in  the  expansion.  On  the  other  hand,  the 
Woodward  synthesis  procedure  is  quite  easy  computationally  but  lacks 
error  estimates  for  points  that  fall  between  the  selected  beam  positions 
for  which  the  pattern  is  exactly  matched.  It  is  thus  impossible  to  pre¬ 
dict  the  required  number  of  array  elements,  or  the  aperture  width,  to 
meet  given  error  specifications,  and  in  fact,  the  question  arises  whether 
the  error  decreases  at  all  for  a  larger  number  of  terms  in  the  approx¬ 
imating  sum. 

Most  of  these  problems  can  be  solved  if  one  uses  trigonometric 
interpolation.  The  technique  is  both  straightforward  and  computation¬ 
ally  easy  and  accurate.  The  estimated  maximum  error  is  of  the  same 
order  of  magnitude  as  that  of  the  Fourier  approximation,  but  the  actual 
error  may  sometimes  exceed  the  Fourier  error  and  sometimes  be 
smaller. 


II.  TRIGONOMETRIC  INTERPOLATION 

The  method  of  interpolation  is  rather  well  known,  but  it  appears 
that  it  has  been  overlooked  in  many  applications.  One  reason  is  that 
often  poor  results  have  been  obtained  from  this  method;  another  is  the 
difficulty  of  inverting  a  large  matrix  and  the  inherent  numerical  inac¬ 
curacies  incurred  in  such  inversion  even  with  the  aid  of  a  computer. 

The  cause  of  both  problems  lies  in  the  fact  that  the  various  terms  of  the 
polynomials  usually  used  for  interpolation,  namely  powers  of  x,  are 
not  mutually  orthogonal.  The  large  errors  of  approximation  can  be  re¬ 
duced  by  introducing  a  non-uniform  spacing  between  interpolation  points 
which  if  properly  chosen  will  result  in  a  Tchebyscheff  polynomial.  But 
the  inversion  difficulty  still  remains. 

In  the  case  of  a  trigonometric  or  exponential  polynomial  the  various 
terms  are  mutually  orthogonal,  not  only  with  respect  to  the  integral  over 
the  period  of  orthogonality  but  also  with  respect  to  the  sum.  Thus  if  the 
functions  and  are  mutually  orthogonal  we  have 

n 

/  d.(x  )  )  =  0  for  i?ik  (1) 

/I _ ,  1  aka  '  ' 

a  -  1 
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where  Xg^  are  the  interpolation  points  covering  the  orthogonality  period  as 
in  the  integration  case.  The  expansion  coefficient  is  given  by 


n 


c.  = 


\  y  (x  ) 


i  ~ 


n 


Z!  “i  '*■>’ 


a=  1 


(2) 


where  is  the  value  of  the  approximated  function  at  x^ .  That  this  is 
the  correct  value  can  be  seen  from  the  following.  Let 


z 

;  =  n 


a  =  1,  2  .  .  ,  n 


(3) 


which  defines  the  interpolation  values  and  their  location.  Multiplying 
both  aides  by  and  summing  over  a,  we  obtain 


N 


j  =  o 


(4) 


a=  1 


0=1 


Interchanging  the  order  of  summation  on  the  left,  namely,  summing  over 
d  first,  we  obtain 


N  n  ,  n  > 


(5) 


j  =  0  0=1 


a  =  1 


using  the  orthogonality  given  in  equation  (1).  Thus, 


(x  )  =  "V  '  y  (x  ) 
1  /  .  j  '  a  /  ,  ^o  i'  o' 

0=1  0=1 


(2a) 
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which  is  equivalent  to  equation  (2). 


In  the  case  of  complex  the  expression  for  c^  will  be  given  by 


c. 

1 


n 


a  =  1 


i^h) 


•  >  a 

where  6.  ia  the  complex  conjugate  of  The  orthogonality  condition  is 
given  by^ 


n 

z 

1 

% 

n 

The  summation  \ 

2 

is  the  normalization  constant  for  the 

a=  1 

vector  rfj^{x). 

For  an  array  of  2N  +  1  isotropic  elements  spaced  \/2  apart,  the  set  of 

functions  forming  our  orthogonal  function  space  are  ?  where  u  =  sin  0. 

The  expression  representing  the  array  is  given  by  J 

N 

P(u)  =  ^  (6) 

n=  -  N 

and  the  excitation  coefficient  a  is  given  by,  ^ 

n 


N  ,  ,  a 

L 

a=  -  N 


(7) 
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where  the  prime  /  ,  implies  that  the  first  and  last  ordinates 
are  taken  with  half  their  value.  Also,  the  first  and  last  excitation  coef¬ 
ficients  are  applied  with  half  of  their  computed  value.  ^  The  reason  for 
this  weighting  factor  which  is  applied  to  the  end  forms  may  be  seen  from 
the  expression 

-i-e- =  sin  nQ  cotO/2  (8) 


which  vanishe 


a  for  9  =  ,  where  m  =  0,  1,  2...  .  If^i  =  ^  then 

n  * 

^  j(i-k)x  J(i-k)^ 


and 


IL,i 


n 


=  0  for  i  k 


(9) 


a  =  -  n 


which  satisfied  the  orthogonality  property  with  the  above  weighting.  One 
actually  may  omit  the  end  element  and  construct  the  interpol|tion  without 
it,  namely,  use  2N  elements  covering  the  interval  of  2it  (  ^  )  of  the 

interpolation  range,  the  reason  being  that  the  value  of  the  ordinate  at  2it 
is  identical  with  its  value  at  zero.  The  sum  corresponding  to  equation 
(8)  will  be 


1-e^ 


for  9 


mrr 

n 


(10) 


Thus 


e  =  0  for  f  k  (11) 


a  =  0 


In  the  case  that  y(u),  the  desired  function,  is  even  or  odd,  the  cor¬ 
responding  expressions  for  the  coefficients  will  be  simplified  and  given  by 


N 


a  =  0 


cos  na 


N 


2 

a  = 
n  N 


(12) 


for  an  even  y(u),  where 


1  N 

y{u)~-r-  a  +a,  cosiru  +  a,  cos  2itu  +  .  .  .  a^.  ,  cos(N  -  l)7ru  +  — ^cos  Nttu 

^''2ol  2  N-1''  2 

(13) 

For  an  odd  y(u) 

N  -  1 

a  =  1 

where 

y(u)  ^  bj  sin  ttu  +  ain  2ff  u  +  ....  +bj^  ^  sin(N  -  1)  iru  (15) 

Note  that  the  cosine  expansion  has  the  same  characteristic  as  the  ex¬ 
ponential,  in  that  the  first  and  last  terms  appear  with  half  value  both  in 
the  expansion  and  in  the  computation  of  the  coefficients  of  excitation  (see 
Appendix  1).  Thus  it  is  seen  that  the  solution  of  the  interpolation  problem 
does  not  involve  any  matrix  inversion,  therefore  ensuring  that  no  loss  of 
accuracy  occurs  throughout  the  computational  process. 


III,  ERROR  ANALYSIS 

The  problem  of  error  estimates  is  becoming  more  important  because 
of  the  stringent  requirements  of  modern  applications,  And  in  view  of  the 
increasing  size  of  arrays  used  the  questions  of  convergence  of  the  represent¬ 
ing  series  are  also  becoming  relevant.  It  was  therefore  considered  of  im¬ 
portance  to  discuss  these  characteristics  as  they  apply  to  the  three  methods 
mentioned  above,  namely,  the  Fourier  approximation,  the  trigonometric 
interpolation,  and  the  Woodward  method,  Since  the  Woodward  method  is 
presently  in  wide  use,  a  particular  attempt  was  made  to  investigate  the  con¬ 
vergence  properties  of  this  method  and  also  to  estimate  its  approximation 
error  as  compared  with  the  Fourier  method. 

The  Fourier  approximation  and  the  trigonometric  interpolation  have 
the  same  convergence  properties  which  can  be  expressed  as  follows:  ^ 

If  f(x)  is  continuous  with  modulus  of  continuity  w(6)  (which  means  that  , 
w(  6  )  =  If  (x,)  -  f  (x.)  fo*"  I  *2  ■  *1 1  ~  ^  ^  ®n  Fourier  n  “ 

partial  sum,  then 

f(x)-s  (x)  <Aw(-^)logn  (16) 

n  n 
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fVi 

If  «(6)  log  6 -♦0  as  6 -'■0,  then  a  (x) -*£  (x)  uniformly.  If  f(x)  has  a  p 
derivative  f^  (x)  such  that  ” 


f**  (x^)  -  f^(Xj) 


:S  K  (x^  -  x^) 


(17) 


for  all  Xj^  and  x^ ,  X.  being  constant,  then 


or 


f(x)  -  flj^(x) 


< 


X.  log  n 

P  +"1 


f(x)  -  s  (x) 

I  “ 


(18) 

(19) 


where  A  is  an  absolute  constant  depending  on  p  alone.  The  same  estimate 
holds  fora  trigonometric  interpolating  sum  of  the  same  number  of  terms. 
Thus  the  rnnvergenr.e  properties  of  the  two  approximation  methods  are 
identical..  This  similarity  holds  as  well  for  piecewise  continuous  function. 


Now  let  us  consider  the  Woodward  approximation.  For  a  discrete 
array  of  2N  +  1  elements  spaced  X./2  the  Woodward  expression  would  be 


N 


P(u)  = 


I 

n=  -N 


sin 


12N_^  . 


£n 


2N  +  1 


n 


sih-^  (u 


2n 

nrrr 
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(20) 


which  when  expanded  is  equal  to 


N 


P(u)  = 


z 

n=  0 


a  cos  niru  +  b  sin  ntru 
n  n 


(21) 


tH 

which  is  a  trigonometric  sum  of  the  n  order.  This  sum  equals  the 
desired  function  at  exactly  2N  +  1  points  spaced  at  distances  of  ^ 

from  each  other.  This,  however,  is  by  definition  an  equispaced  trigono¬ 
metric  interpolation  and  as  such  should  have  again  the  same  covergence 
properties  as  the  Fourier  expansion.  Thus  the  convergence  properties 
of  all  three  methods  are  the  same. 

The  next  question  is  that  of  the  error  estimate.  This  problem  is 
very  difficult  if  a  useful  estimate  is  wanted,  namely,  within  a  few  percent 
of  the  actual  error.  This  is  the  reason  why  equations  (16)-(19)  cannot 


I 
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serve  this  purpose  and  are  good  only  as  a  convergence  criterion  when  a\ 
large  number  of  terms  are  involved.  The  difficulty  lies  in  finding  the  lo-W^st 
upper  bound  on  the  constant  A  or  A  in  each  of  the  cases  mentioned  above; 

A  constant  which  is  close  to  this  oprimum  is  available  for  the  Fourier  ap- 
proximation  of  continuous  functions  and  is  given  by^  \ 

\ 

A=-\  (22) 

TT 


SO  that  the  error  estimate  in  the  maximum  absolute  deviation  sense  is 
given  by 


E 


log  n 


(23) 


where  n  is  the  order  of  the  corresponding  Fourier  partial  sum,  and  Mp  is 
the  maximum  value  of  the  p*^  derivative  of  the  approximated  function.  This 
estimate  seems  to  hold  even  when  applied  to  functions  containing  discontinuities 
except  in  the  vicinity  of  the  discontinuity  where  the  occurrence  of  the  Cibb'e 
phenomenon  will  produce  a  deviation  of  approximately  9  percent  of  the  total 
jump  in  the  function.  This  deviation  is  independent  of  n .  This  phenomenon 
can  be  eliminated  by  the  modification  of  the  Fourier  coefficients  in  a  manner 
which  is  equivalent  to  a  Fejer  summation,  ’ 

For  the  trigonometric  interpolation  process  such  minimal  constant  is 
not  available,  but  it  can  be  shown  that  the  error  in  the  interpolation  method 
is  very  close  to  that  of  the  Fourier  approximation.  It  will  be  shown  that 
the  difference  between  the  two  is  in  fact  that  the  Fourier  partial  sum  of  n^“ 
order  completely  neglects  higher  order  components,  while  in  the  interpola¬ 
tion  method  the  higher  order  components  get  reflected  and  add  on  to  the 
lower  order  components.  ^  However,  for  a  good  Fourier  approximation 
the  deviation  of  the  n^“  partial  sum  from  the  desired  function  la  very  small; 
this  implies  that  the  higher  order  components  are  negligible,  and  hence  it 
makes  little  difference  whether  their  amplitudes  are  recorded  or  not.  Also, 
one  should  remember  that  although  any  modifications  of  the  Fourier 


Reference  4,  page  105. 
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components  increase  the  mean  square  error,  they  do  not  necessarily 
increase  the  maximum  deviation;  in  fact,  in  some  cases  they  decrease 
it.  One  such  example  is  the  tapering  of  higher  order  Fourier  components 
which  was  shown  to  lead  to  a  srnaller  maximum  deviation  than  that  of  the 
equivalent  Fourier  partial  sum.  Let  us  consider  in  detail  the  effect  of 
substituting  the  interpolation  method  for  the  Fourier  approximation.  If 
the  approximated  function  f(x)  has  only  the  first  N  Fourier  components, 
both  methods  will  yield  the  function  exactly.  In  the  case  of  higher  order 
components  they  will  be  recorded  in  the  following  way:  consider  the 
(N  +  k)^“  Fourier  component  of  the  m*"  term  in  the  interpolation  sum 
and  assume  we  have  an  odd  function  such  that  only  sines  are 

present. 

^  ^  =  sin  (N  +  kjxsin  mx  (24) 


=  (sin  Nx  cos  kx  +  cos  Nx  sin  kx)  sin  mx 


for 


X  =  a 


IT 


tt  =  1,  2  ...  N  -  1 


j  N  +  k  ,  .  ,  air  .  . 

^m  ~  {coa  air  am  k  — nix 


til  til 

Now  consider  the  (N  -  k)  Fourier  cpmponent  of  the  term 


(25) 


sin  (N-k)  X  sin  mx 


(26) 


The  mean  square  error  is  equal  to 


a' 


n 


where  aj^  is  the  Fourier  coefficient  and  =  a  +  e  ,  and  f(x)  is  the 
approximated  function,  Substituting  for  we  ^et  ^ 


1  I 

"ZiT  I 


f(x) 


'  dx  + , 


(a  +  «_)(€_  -  a  ) 
n  n  n  n 


where  the  equality  sign  holds  only  for  =  0  • 
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results  in 


which  for  x  = 


I  N  -  k  ,  .  ,  Off .  .  ,  .  j  N  +  k 

9  =  -(cos  air  sin  k  sin  (mx)  =  -o 

m  '  N  '  '  '  m 

Thus  it  is  seen  that  the  higher  odd  harmonic  components  of  the  function 
f  (x)  will  be  recorded  with  reversed  sign  as  lower  components.  The  even 
harmonics  are  recorded  with  the  same  sign  as  can  be  seen  from  the  fol¬ 
lowing  equation, 

N  +  k  air  N  -  k 

=  cos  (N  +  k)xco8  mx  =  cos  air  cos  k  =  coa(N  -  k)x  cos  mx  =  W 
m  N  m 

for  X  =  (Z7) 

As  mentioned  above,  if  the  higher  order  components  are  negligible  in 
magnitude,  their  effect  will  be  just  as  negligible  in  the  Fourier  series  as 
it  is  in  the  interpolation  case.  In  the  case  that  these  components  are  sig¬ 
nificant  the  implication  is  that  the  Fourier  approximation  will  have  a 
significant  error.  In  this  case  the  interpolation  will  also  suffer  from  a 
similar  although  not  identical  error,  and  the  maximum  deviation  error  will 
be  somewhat  smaller  in  some  cases  and  somewhat  larger  in  others,  but 
it  will  be  of  the  same  order  of  magnitude  as  that  of  the  Fourier  approxima¬ 
tion. 


Let  us  now  investigate  the  error  involved  in  the  Woodward  technique. 
As  mentioned,  the  Woodward  technique  is  an  interpolation  method  with  its 
interpolation  points  slightly  closer  together  than  in  the  method  previously 
discussed.  Thus  for  an  N*”  order  trigonometric  sum  the  spacing  la 

.  in  the  Woodward  as  compared  to for  the  former  one.  One  can 
Intuitively  say  that  the  two  methods  shoiad  yield  almost  the  same  results 
when  N  is  relatively  large,  and  assuming  that  f(x)  is  continuous.  We  shall 
try  to  show  in  a  more  rigorous  fashion  that  the  two  interpolations  are  quite 
similar  and  also  indicate  in  what  way  they  differ. 

To  simplify  the  discussion  we  will  assume  that  f(x)  is  an  even  function 
of  x;  this  does  not  result  in  loss  of  generality  since  a  similar  proof  may  be 
used  for  an  odd  function  and  a  sum  of  the  two  can  represent  any  arbitrary 
function.  Since  f(x)  is  even,  it  is  represented  by  a  summation  of  cosines 
alone.  Thus,  let  the  regular  interpolation  sum  discussed  previously  be 
represented  by  Pj(x)  and  the  Woodward  sum  by  hence 

N 

p.w  =  ^  cos  nx  (28) 

n=“'o 
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where  Pj(x}  =  f(x)  at  x  =  oi  =  0,  +  1,  .  .  .  i  N 

N  ■ 

P2{x)  = 

n  =  0 

where  =  f(x)  at  x  =  a  ^  a  =  0,  i  1,  .  . .  +  N 


Zb  COB  nx 


(29) 


Let  us  now  substitute  powers  of  cos  x  for  cos  nx  such  that 

N  N 


P^(x)  =  'y  '  cos  nx  a  ^  |  cos^ 


and 


0 

N 


0 

N 


Pp(x)  =  y  '  b^  cos  nx  =  y  cos*'  x 


(30) 


(31) 


Now  substituting  y  for  cos  x  we  obtain, 

N 

0 

and 

0 

where 


(32) 


(33) 


0  <  X  £  TT 


+  1  ay  a  -  1 


We  now  have  two  polynomials  of  y  interpolating  the  function  f(y)  where 
y  =  cos  irx.  This  amounts  to  a  Lagrangian  interpolation  whose  char¬ 
acteristics  are  well  known  and  where  the  error,  or  remainder,  is  given 
by 


R(y)  = 


-idixL. 

(N  +  1) 


^N+1^^) 


(34) 
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where 


(35) 


w(y)  =  (y  -  Pq)  (y  -  p^)  . .  .  (y  -  p^) 

and  where  the  p^j  are  the  interpolation  points.  At  these  points  the  function 
is  equal  to  tl.ie  polynomial  and  therefore  the  error  must  vanish.  The  value 
4  represents  a  point  close  to  the  middle  of  the  interpolation  interval  of 
P(y),  and  is  connected  with  the  location  of  p^^- 


Now  let  us  cornparethe  two  remainder  expressions 

n  ,  .  (y-l)(y-cos^)(y-cos  -^)  .  .  ■  (yi-1) 

Ri(y)-  (N  +  i)  :  ^  ^^1^ 


and 


2it 


2irN 


R,(y)  = 


(y-l)(y-cos^^^)...  (y-cos-^^^) 
(N+  1)  1  ^ 


(36) 


(37) 


The  difference  between  the  expression  is  in  4  and  in  the  p  .  For  a 

N+1 

fairly  smooth  function  the  variations  of  f  (x)  are  limited.  The  change 
in  4  is  directly  related  to  the  distribution  of  p^^ ,  and  a  small  shift  in 
those  points  will  produce  a  correspondingly  small  change  in  4  •  E<xag> 
gerating  this  shift  by  assuming  that  all  points  moved  by  an  amount: 


TT  ZlT  _ IT 

N  ■  In  +  1  ""  N(2N  +  1) 


we  obtain  that  the  change  in 


^N+1 


(4)  is  given  by 


(38) 


-  f^'*'^42)|  s 


,N+2 


N(2N+  1) 


(39) 


The  difference  in  the  p  's  is  reflected  in  the  values  of  w(y)  .  Now  for 
most  of  the  interval  of  Approximation  the  zeros  of  are  more  closely 

spaced  than  that  of  W2(y);  hence,  the  magnitude  of  oo2(y)  S  wi(y)  for  this 
range.  However  the  error  will  increase  in  the  \acinity  of  y=  -1  (x-  ir), 
since  a  zero  located  there  was  shifted  to  x  =  ^^+1  '  value  of  the 

error  may  be  calculated  if  desired  by  computing  the  total  effects  of  the 
zero  shifts.  From  a  practical  point  of  view  this  often  is  of  no  importance 
since  the  element  factor  may  produce  a  zero  at  precisely  that  location, 
which  is  the  end  of  the  visible  range.  Thus  for  the  range  of  interest  and 
for  N»1 
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We  can  therefore  see  that  for  relatively  long  arrays  the  error  in  a  Wood¬ 
ward  synthesis  is  not  substantially  different  from  the  Fourier,  except  at 
the  end  of  the  visible  range.  If  the  array  is  not  subject  to  scanning,  this 
may  not  be  objectionable  at  all. 

For  continuous  apertures  the  approximation  characteristics  of  the 
Woodward  method  are  of  similar  nature,  and  it  can  be  shown  that  the 
approximating  sum  converges  uniformly  to  the  desired  (continuous)  pat¬ 
tern.  (For  a  discussion  on  the  continuous  aperture  approximation  see 
Appendix  11.  ) 


IV.  EXAMPLES 

In  order  to  provide  a  comparative  example  for  the  above  methods  of 
approximation  a  simple  pattern  was  chosen  which  is  easy  to  obtain  by  means 
of  any  of  the  above  methods,  namely  a  sector  beam  of  ±30°  beamwidth. 
Since  the  beam  contains  discontinuities,  it  is  a  good  example  of  the  Gibb's 
phenomenon  as  well.  If  we  compare  Figure  1  with  Figure  Z,  we  can  see 
that  the  Interpolation  method  yields  excellent  results  except  for  the  over¬ 
shoot  area  which  seems  to  suffer  even  more  than  the  original  Fourier  ap¬ 
proximation.  This  tends  to  confirm  the  original  contention  that  in  some 
cases  errors  may  be  worse  than  those  of  the  Fourier  method  and  in  some 
cases  better.  Thus  it  seems  that  in  the  continuous  range  the  Interpolation 
method  results  in  a  smoother  pattern  and  in  the  discontinuous  it  is  more 
oscillatory.  Turning  now  to  the  Woodward  approximation  (see  Figure  3) 
it  is  seen  that  the  error  is  substantially  larger  than  in  the  former  cases 
at  the  end  of  the  visible  range  while  it  is  quite  similar  at  the  origin,  in 
fact  as  the  remainder  expression  Indicated  the  error  at  the  origin  is  smaller 
than  either  the  error  of  the  Fourier  method  or  of  the  interpolation  techniques. 


INTERPOLATION 
EXACT  - 


V. 


CONCLUSIONS 


A  method  has  been  presented  by  means  of  which  an  approximation  to 
any  desired  pattern  can  be  synthesized  having  error  characteristics  very 
similar  to  those  of  a  Fourier  approximation.  The  method  is  simpler  com¬ 
putationally  than  the  Fourier  and  consequently  is  subject  to  fewer  computa¬ 
tional  errors  as  well  as  fewer  expenses. 

The  convergence  and  error  characteristics  of  the  Fourier,  the  Wood¬ 
ward,  and  the  interpolation  methods  have  been  presented  and  an  example 
has  been  computed  to  demonstrate  the  error  behavior.  The  results  indicate 
that  all  three  methods  are  adequate  in  many  instances  but  in  some  the 
Woodward  method  may  not  be  satisfactory. 


APPENDIX  1 

Computation  of  the  Coefficients  of  a  Fourier  Interpolation  Sum 
of  6  Terms  Approximating  a  Sector  Beam  (see  Figure  2). 

The  first  coefficient  is  given  by 

5  . 

So  =  2  ^  ^  ^  2  1  +  1)  =  0,  5 

O' 

The  coefficients  of  cos  niru  are 

a^  -  '  y^cos  (0.  5  +  cos  0.  2ir  +  cos  0.  4^ )  =  0.  646 

32  =  '  y^  cos  =  "I"  (0.  5  +  cos  0,  4it  +  cos  0.  Sir )  =  0 

a^  =  '  y^cos  =  ”1“  (0.  5  +  cos  0.  bir  +  cos  1 .  Zir )  =  -0.  247 

a^  =  ^  y^  cos  ~“f'  (0*5  +  cos  0.  8ir  +  cos  1.6it)  =  0 

Sg  =  ~Y  x-|-^  I  y^cos  =  -|-  x-|'  (0.  5  +  cos  IT  +  cos  2tt  )  =  0.  100 

Note  that  a^  and  are  divided  by  2.  The  resultant  expression  for  the  poly¬ 
nomial  is  given  by 

F(u)  =  0.  5  +  0.  646  cos  21X0  -  0,  247  cos  .3tru  +  0.  100  cos  5tu. 
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APPENDIX  II 


Approximate  Characteristics  of  Woodward  Synthesis  for  Continuous  Apertures. 


The  Woodward  synthesis  is  closely  related  to  the  sampling  theorems 
in  networks,  and  certain  recent  mathematical  results^  obtained  in  that  field 
can  be  applied  when  properly  interpreted.  For  an  infinite  aperture  the 
pattern  can  be  obtained  exactly  as  long  as  the  transform  pair  exists. 


00  00 

G(u)  =  I"  F(x)e-**“dx  ;  F(x)  =  J  G(u)  e“'^'*'^du  (II- 1) 

-  00  -00 

where  u  s  tt  sin  0  and  x  is  measured  in  half  wavelength  units  (-j-^  where 
d  has  the  dimension  of  length),  F(x)  represents  the  excitation  while  G(u) 
-vepresents  the  far  field  pattern.  The  case  of  practical  interest,  however, 
involves  apertures  of  finite  size.  One  can  obtain  an  upper  bound  on  the 
error  of  the  Woodward  approximation  for  such  cases,  and  one  can  also 
show  that  this  error  goes  uniformly  to  zero  as  the  aperture  Increases  in¬ 
definitely.  The  error  depends  on  the  behavior  of  the  part  of  the  excitation 
distribution  that  la  omitted  due  to  the  finiteness  of  the  available  aperture. 

We  will  assume  that  this  excitation  attenuates  at  a  certain  rate  (if  it  did  not 
attenuate  the  transform  pair  of  equation  (II- 1)  would  not  exist)  given  by 


F(x) 

F'(x) 


<M  j  x 


-k 


X 


x>W 


(11-2) 


where  F'(x)  is  the  derivative  of  F(x),  W\  is  the  aperture  width,  and 
M2>0,  M2\>0,  k^  >  2,  k2  >  2  .  Another  constant  representing  the 
characteristics  of  the  pattern  both  in  the  visible  and  invisible  ranges  is 
given  by 


^3=  jk 


00 


du 


(II-3) 


The  behavior  of  G(u)  outside  the  visible  range  may  be  arbitrarily  specified, 
and  since  it  represents  the  reactive  energy,  it  may  be  advantageous  to  re¬ 
quire  that  it  be  small.  One  must  remember  though,  that  its  specification 
will  affect  the  conditions  given  by  equation  (II-2)  which  affect  the  error  ex¬ 
pressions.  The  error  bound  is  given  by 
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N  -  1 


G(a)-  ^ 


n  =  N  -  1 


1  -  k, 


k,  +  4  2  -  k, 

^ — r-  W  ^ 


+  W' 


5M 


M-  + 


3  8(k2-  1) 


2  -^2'' 
Q-  -  tir  « 


w 


-i2— -  +  -tH-t  (2Nir  -  1) 
2Nit  -  1  Nit  -  1  J 


(11-4) 


From  equation  (11-4)  one  can  see  that  the  third  term  on  the  right  ha*^  aide 
goea  uniformly  to  zero  when  N  goea  to  infinity,  and  W  remalna  •  " 

ting  N>W  Impliea  that  beams  are  pointed  ’ 

the  presence  of  some  reactive  energy.  However,  this  does  not  imply  *  P  ’ 
gain;  it  juat  requires  that  the  reactive  components  of  the  excitation 
be  properly  adjusted.  Under  such  conditions  the  error  bound  is  given  by 
the  first  two  terms  alone,  namely, 


I-kj  ^2  +  4  2-1^2 

Error26MjW  W 

which  goes  to  zero  as  Wi  or  the  aperture  width  increases  indefinitely. 
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